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The dynamical development of the cooling and hadronizing quark-gluon Plasma (QGP) is studied 
in a simple model assuming critical fluctuations in the QGP to Hadronic Matter (HM) and a first 
order transition in a small finite system. We consider an earlier determined free-energy density curve 
in the neighbourhood of the critical point, with two local minima corresponding to the equilibrium 
hadronic and QGP configurations. In this approach the divergence at e = eliminates fluctuations 
with negative or zero energy. The barrier between the equilibrium states is obtained from an 
estimated value of the surface tension between the two phases. We obtain a characteristic behavior 
for the skewness and the kurtosis of energy density fluctuations, which can be studied via a beam 
energy scan program. 
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I. INTRODUCTION 

In central heavy ion collisions fluctuations may occur 
due to critical phenomena arising from a phase transi- 
tion in the Equation of State (EoS). Fluctuations arising 
from the initial asymmetry are smaller in these head-on 
reactions. In the neighborhood of the critical point the 
shear viscosity of the quark-gluon Plasma (QGP) is be- 
coming small [l[ , which also facilitates the appearance of 
fluctuations. 

The flow effects depend on the initial state profile and 
on the viscosity. Turbulence appears only for small vis- 
cosity, which indicates the critical point of the matter [l[ , 
and it is a sensitive measure of viscosity and its minimum 
at the critical point. 

Molecular dynamics and Fluid Dynamics simulations 
of heavy ion collisions suggest that collective flow asym- 
metries can be measured Q, both if the global asymme- 
try or random flow arising from the initial state trans- 
verse momentum fluctuations or longitudinal center of 
mass rapidity, ycM, fluctuations are causing it. These 
alternative sources may also lead to specific statistical 
characteristics as discussed in ref. [|| . These phenomena 
contribute to a spatial spread of the matter and energy 
density variations are present even if we do not have a 
phase transition in our EoS. 

Here, on the other hand, we study random fluctuations 
of thermodynamical origin caused by the phase transition 
based on the considerations described in ref. [J- The 
field is under intensive recent study, and several aspects 
of the phase transitions have discussed the non-Gaussian 
fluctuations 

For a realistic reaction model we have to describe the 
final stage of the reaction also. In order to study the 
properties of the phase transition we need to have suffi- 
cient energy to form QGP in a sufficiently large volume, 
and then the system must hadronizc. At relatively low 
energies this hadronization might happen well before the 
freeze out (FO) stage, and then we have only rather in- 



direct information about the phase transition. 

In our recent studies Q we use a Multi Module Model 
or Hybrid Model approach to describe high energy heavy 
ion collisions in the RHIC and LHC energy range. Then 
from the locally equilibrated QGP we have to form 
hadrons. We do not assume that the hadronization 
happens in chemical equilibrium as this would take too 
long time @ and would not allow for baryons of high 
strangeness. 

In high energy heavy ion reactions around 50 
GeV/nuclcon or above when QGP is formed, we have 
low mass quarks (5-10 McV) and massless gluons in 
large numbers and the conserved baryon charge has a 
minor effect. On the hadronic side also mainly mesons 
and baryon-antibaryon pairs are formed. The observ- 
able phase transition happens when the plasma expands, 
hadronizes, and freezes out. At this stage the system may 
be in the vicinity of the critical point so critical fluctua- 
tions may appear. In a finite volume this would show up 
as energy density fluctuation, which would then lead to 
charged hadron number fluctuation also (and much less 
in net baryon charge fluctuations). 

In the present simple model we assume two coexisting 
phases in a finite system near a first order phase transi- 
tion. We study the abundance of the energy density dis- 
tribution of the two phases in the mixed phase domain 
in terms of the volume ratio. 

This dynamical development is not directly observ- 
able, nevertheless, we can observe the FO moment, which 
should be similar for different locations in the system. 

At the same time we can vary the beam energy, which 
shifts the FO point versus the critical point, so a beam 
energy scan can provide us the series of information we 
arc interested in. 
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II. CRITICAL FLUCTUATIONS 

Following ref. Q we briefly present the way to describe 
critical fluctuations following the ideas of the Landau- 
Ginsburg theory for critical phenomena. 

In case of QGP to hadronic matter (HM) transition 
the essential difference is that in phase equilibrium the 
energy density of the HM phase is much lower than that 
of the QGP phase. As a consequence the usual 4th or- 
der polynomial to describe the free energy of the system 
is not realistic as it would lead to considerable popula- 
tion of negative energy density states, which would be 
unphysical. 

This problem was solved in ref. where the 4th 
order polynomial approach was substituted by a Laurent 
series. 

Using the simple bag model EoS the equilibrium value 
of the energy density in the low-temperature, low-energy- 
density phase (HM) denoted by e^, and the equilib- 
rium value of the energy- density in the high-temperature, 
high-encrgy-density, QGP, phase denoted by e q can be 
calculated as 
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To study the fluctuations we find the free energy, /(e), for 
arbitrary values of e, and not just for the energy densities 
e q and eh- For these equilibrium points: 



f(e q (T)) = - Pq (T) 



-(37T 4 -34T 4 ), (3) 



f(e h (T)) = - Ph (T) = - 
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Following the Landau theory we approximate now the 
free energy density as a polynomial in the neighborhood 
of an eo energy density (eo € [e^, e 9 ]), where it has a local 
maximum. In order to obtain the required divergence at 
e = 0, a slightly modified functional form is assumed: 

/(e) = h + — + K 2 {e - e ) + K 3 (e ~ e ) 2 + K±(e - e ) 3 . 
e 

(5) 

The constants, f±, K\, K 2 , K 3 , K± can be determined 
from thermodynamic considerations as done in ref. 
One will obtain 
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K-2 = (7) 



K 3 = - 
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In Eq. ((6]) a represents the surface tension of a hadronic 
bubble and £0 is the characteristic size of a hadronic 
droplet. Once the Ki values are expressed as a function of 
e, eo, e q , eh, £o> a values, the unknown /1 and eo param- 
eters are obtainable at any temperatures from Eqs (|3]) 
and (U). Following our previous work Q, we have used 
the T c = 0.169MeV, £ = 3/m and a = 0.05GeV/fm 2 
values in all our calculations. 

For temperatures in the vicinity of the critical temper- 
ature the free energy density curve as a function of the 
energy density is illustrated in Fig. [TJ 
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FIG. 1: (color online) The free energy density as a function 
of energy density e, for T = 0.95T C (dashed), T — T c (full 
line) and T = 1.05T C (dot-dashed). 



Once the free energy curve is known, one can esti- 
mate the probability density of finding the system in a 
state with energy density e: P(e) oc exp(— /3F(e)), where 
F(e) = ttf(e), with f2 the volume of the created QGP. 
On Fig. [5] we plot at the critical temperature the char- 
acteristic P(e) / P(e q ) curves, considering different vol- 
umes for the QGP (Q. = 10 fm 3 , 50 fm 3 , and 500 fm 3 
values). Also, in Fig. [3] we show the P(e)/P(e q ) curves 
for fl = 500 fm 3 and the same temperatures as in Fig. [T] 

It is also important to mention that different thermo- 
dynamical parameters (especially intensives and exten- 
sive ones) do not have to show the same critical fluctu- 
ation properties, so we have to study the fluctuations of 
several parameters. Furthermore, the statistical physics 
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estimates assume a single thermal source at or near the 
critical point, while we estimate here also the effects of 
spatial fluctuations, which arise from a dynamically ex- 
panding fluid flow even in the least fluctuating configu- 
ration. 
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FIG. 2: (color online) The relative probability of finding a 
state of a given energy density e for T = T c , in a system 
of volume Q = 10 fm 3 (dashed), Q = 50 fm 3 (full line) and 
n = 500 fm 3 (dot-dashed). 



We assume central collisions only, to avoid the effects 
from azimuthal flow asymmetries and from particle emis- 
sion fromprojectile and target residues (spectator evap- 
oration) Q. 

In a theoretical approach we can assume a spatial dis- 
tribution of the x quantity. For the variable x the av- 
erages and various order moments distributions can be 
written as 

(x n ) = [ x n P(x)dx, (12) 



MM = ((x-(x)r)= (13) 

where P{x) is the spatial distribution weighted, e.g., by 
the baryon charge density in the center of mass frame 
(CF). The spatial variance, the skewness and the kurtosis 
can be obtained from these moments: 

Ax= ((x- (x)) 2 ) = M {2) , (14) 

(Ax) 3 /2 (M( 2 ))3/2 ' {L0) 
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FIG. 3: (color online) The relative probability of finding a 
state of a given energy density e for the T = 0.95T C (dashed), 
T — T c (continuous) and T = 1.05T C (dot-dashed) tempera- 
tures. The of volume of QGP is Q = 500 fm 3 . 



III. HADRONIZATION AND EXPANSION 
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By using these averages, first we can calculate spe- 
cific extensives, which are governed by strict conservation 
laws. The total baryon charge, energy and momentum 
conservations are governed by the continuity equation 
and by the relativistic Euler equation. 
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and as a consequence, the total momentum in the CF 
should remain zero during the development, while the 
average specific energy per net nucleon number 



,CF\ 



llll 



T 



const. 



(19) 



should remain constant in CF. 

The total baryon number, N to t, is exactly conserved 
in the reaction. At the same time the average baryon 
charge density is decreasing. 



The dynamically developing flow pattern leads to a 
spatial distribution of all thermodynamical quantities, 
while the system expands rapidly. Finally the super- 
cooled QGP can hadronize rapidly and almost simulta- 
neously it freezes out. This final stage of the reaction can 
be described by a non-equilibrium model. 



IV. SKEWNESS AND KURTOSIS 

In this section we study the skewness and kurtosis of 
the specific energy density and consequently the charged 
particle densities according to Eq. (fT5j) and Eq. (fT6|) . We 
determine these quantities as a function of the systems 
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temperature and also as a function of the volume abun- 
dance (rh) of the hadronic matter. We assume that in a 
rapid transition, where critical fluctuations dominate and 
the two phases are not separated these two phases are in 
thermal equilibrium. The temperature decreases rapidly 
starting from QGP where T > T c until the hadronization 
completes at T < T c . As the phases are not separated 
the simplest estimate for any temperature T, is that the 
volume abundance of the hadronic matter is: 



P(e h ) 



P(e q )+P(e h y 



(20) 



where and e q are the energy densities of the pure 
phases defined in eqs. (TTJ [5]), and the probability den- 
sities, P(e), are defined in section [TTJ This relation makes 
a one to one correspondence between the volume abun- 
dance and the equilibrium temperature in a rapidly ex- 
panding and cooling system during the process of a phase 
transition. 

The skewness, Fig. 01 is first negative (indicating a 
longer tail on the lower energy side), then at 80 % HM 
(rh = 0.8) volume abundance it turns into positive (in- 
dicating a longer tail on the high energy side). The 
hadronization can be parametrized both function 
of the volume abundance of the growing hadronic phase 
or the decreasing equilibrium temperature of the system 
with critical fluctuations. 

In Fig. [51 we can see that the kurtosis is positive at 
first, then turns to be negative (the distribution becomes 
wider) in the phase transition domain, while it becomes 
positive again as the phase transition completes. The 
minimum of kurtosis is at 80 % HM volume abundance. 
Notice that the kurtosis is increasing much sharper on 
the hadronic side, which is a clear consequence of the en- 
ergy difference between the two phases. This asymmetry 
appears as a result of the Laurent series expansion, and it 
would not show up with the usual 4th order polynomial 
approximation! 

In a multi-module or hybrid-model construction (e.g 
where the PIC hydro stage is attached to a parton and 
hadron cascade model PACIAE), the flow features are 
matched [l5| to a subsequent dynamical model, which de- 
scribes dynamical, non-equilibrium, rapid hadronization. 
This type of models can describe realistically the statis- 
tical properties in a dynamical phase transition, which 
determine the hadron distribution in the final stage. 

This stage would then explicitly describe the random 
fluctuations arising both from the phase transition and 
the flow dynamics. 

In ref. fl6j a mixed particle method is introduced, 
which could separate the fluctuations arising from local 
critical fluctuations. The "mixed events" are actually 
eliminating two particle correlations, and only the sin- 
gle particle distributions remain. Thus for central events 
these are mainly local correlations, which may arise from 
local fluctuations caused by energy and baryon charge 
clustering in a phase transition. The method separates 
the consequences of such correlations. 
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FIG. 4: (color online) Skewness as a function of the volume 
abundance of the hadronic matter (denoted as r^, where 1 
represents complete hadronization). The temperature scale is 
also indicated for clarity, the identifiers represent increments 
of 0.1 MeV in T. Results for ft = 500/m 3 . 
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FIG. 5: (color online) Kurtosis as a function of the volume 
abundance of the hadronic matter (denoted as rh, where 1 
represents complete hadronization). The temperature scale is 
also indicated for clarity, the identifiers represent increments 
of 0.1 MeV in T. Results for Q = 500/?n 3 . 



This method can be used both in hybrid model calcu- 
lations and in experiments, to separate the fluctuation 
effects form the collective flow and the phase transition 
dynamics. 



V. CONCLUSIONS 

In this work we reiterated earlier results on critical 
fluctuations Q in the phase transition between quark- 
gluon plasma and hadronic matter. This model is spe- 
cific because of the large difference of the energy density 
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of the two phases, where e q ^> e%. Here this work is 
extended to the evaluation of frequently used statisti- 
cal parameters like the Skewness and Kurtosis of typical 
parameters like charged particle multiplicities. In con- 
tinuum models the energy density fluctuations or baryon 
charge fluctuations could carry the same role, however at 
ultra relativistic energies a large number of antiparticles 
are created therefore the energy density fluctuations are 
better representing the total charge particle multiplicity 
fluctuation. 

Our model is furthermore extended to study the dy- 
namical change of these typical parameters during the 
hadronization process, where the development is very 
characteristic and informative for the phase transition 
we study. 

In experiments one can measure these parameters (to- 
gether with all other measurable quantities) at the Freeze 
Out (FO) time (or at the FO hypersurface) . Luckily the 
recent RHIC Beam Energy Scan program, scans the same 



statistical parameters at a series of different beam ener- 
gies, which also provide a set of different FO points map- 
ping the dynamics of the phase transition or the most 
interesting part of it. The FO point at high energies 
where QGP is formed may be simultaneous to the phase 
transition, as hadronization must always take place in 
these reactions. At lower beam energies when QGP is 
not formed the FO is in the hadronic phase. At very 
high energies where high energy density QGP is formed 
in a large volume, the system may have sufficiently large 
volume and large energy density that FO happens later 
then the hadronization, i.e. in the pure hadronic state. 

The dynamical changes of the statistical parameters 
this way can provide valuable information at what stage 
of the phase transition we are at a given reaction. 
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